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Abstract 

We study the dynamics of one-dimensional discrete models of one-component 
active medium built up of spatially inhomogeneous chains of diffusively coupled 
piecewise linear maps. The nonhomogeneities ("defects") are treated in terms of 
parameter difference in the corresponding maps. Two types of space defects are 
considered: periodic and localized. We elaborate analytic approach to obtain the 
regions with the qualitatively different behaviour in the parametric space. For the 
model with space-periodic nonhomogeneities we found an exact boundary sepa- 
rating the regions of regular and chaotic dynamics. For the chain with a unique 
(localized) defect the numerical estimate is given. The effect of the nonhomegeneity 
on the global dynamics of the system is analyzed. 

Key words: coupled map lattices, nonhomogeneities, chaos, Chebyshev 
polynomials 

PACS: 05.45.Ra, 05.45.Jn, 05.45.Pq, 05.45. Gg 



1 Introduction 



One of the effective methods of analysis of several systems is based on their 
representation by discrete models. For this one can use the space discretization 
as well as the time one. Using the space discretization the initial (continuous) 
system is modelled by a finite set of elements that are coupled according to 
some rules, which follows from the nature of the system of interest. Each 
element of these then occurs to be some (small) dynamical system which 
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may be described by a (small) set of dynamical variables. If the time is also 
discretized. i.e. if the evolution of these dynamical subsystems is considered 
at discrete time-steps, then such subsystems are treated as maps. In this case 
spatio-tempotal models are said to be lattice of coupled maps. 

The examples of application of the space-time discretization in studies of the 
nonlinear media are numerous. Here we mention the problems related to chaos 

and noise in distributed systems [1,2,3,4], to description of the excitable media 
[4,5] and some other similar problems (see e.g. [6,7,8,9,10,11,12,13,14] and 
references therein). Moreover, any numerical analysis of a continuous system 
has in fact to be reduced to the study of space-time discretized model. 

In principle, a variety of different space-time lattices may be used. The mostly 
studied are the models where only the nearest-neighbour (local) coupling is 
taken into account (see e.g. [8,9,10,12,13] and references therein). In other 
models the global coupling, i.e. when all elements of the system are coupled 
with each other, is assumed [15,16,17]. In lattices with the local coupling, as a 
rule, this is considered as a diffusional one [1,2,4,8,9,10,11,12,13]. Mainly, these 
models are used to analyse the spatio-temporal chaos, pattern formation, self- 
organization phenomena, etc. [18,19,20] 

If all the elements are the same, then such lattices are called homogeneous. 
If, however, elements are different, then the lattice is not a homogeneous one, 
and its investigation is much more complex. With applied point of view the 
nonhomegeneity may be two-fold: either elements of the system have differ- 
ent maps, or they have functionally identical maps but with different control 
parameters. The nonhomogeneities are usually said to be defects. The relative 
position of defects may vary from a single localized inclusion to a space- 
periodic ordering. 

In most of studies of coupled map lattices, the case of homogeneous lattices 
(and mainly numerically) has been addressed. However, it is clear, that the 
assumption of the homogeneity is a crude simplification. Therefore it is rather 
important to investigate the sensitivity of the lattice dynamics to certain non- 
homogeneities (see, e.g., [21,31] and references therein). 

The present study is devoted to the spatially inhomogeneous ID lattices (in 
fact, chains) of diffusively coupled piecewise-linear maps, which refer to some 
models of the one-dimensional excitable media [22]. The defects are deter- 
mined as maps with different parameter values. Using the exact calculation 
of the Lyapunov exponents we analyse the behaviour of the lattice with a lo- 
calized defect and with the space-periodic nonhomegeneity. We also describe 
the structure of the phase space of such systems. 
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2 Models of spatially inhomogeneous excitable media 



Here we consider the dynamics of two models of inhomogeneous active media: 
one-dimensional (annular) chain of coupled maps with periodic nonhome- 
geneity and with one localized defect. To introduce the notations, we remind 
that the term "map" is implied transformation of some interval / on itself: 
Ta : / — >■ /, where / = G R"*^, x i— > G(a;, a), x G I, G{x,a) is some 

function and a is a parameter. In terms of iterations k this may be written as 
x{k + 1) = G{x{k),a), x E I. The nonhomegeneity means nonuniformity of 
the parameter a along the chain of maps. The functions G are chosen in such 
a way that they model the one-component active medium. 

2.1 Homogeneous medium 

Consider the system of N diffusively coupled maps 

Xn{k + 1) G{xn{k), a, 7) + e(xn+i{k) - 2xn{k) + Xn-l{k)) (1) 



with periodic boundary conditions, Xn^N{k)=Xn{k), where n = 1,2, ...,N 
being the discrete space coordinate. A; = 0, 1, 2, ... is the discrete time, e > 
is the (constant) diffusion coefficient, and a > 0, 7 are the parameters of the 
map. This choice of the coupling means that the spatial and temporal actions 
in map (1) occur simultaneously. Let G be the piecewise linear function: 



(1 — 2a)x — 7, 



G(x, a, 7) = a; + aF{x) — 7 = < 



X < 1/2; 



^ (1 - 2q;)x - 7 + 2q;, x > 1/2. 



(2) 



This function G{x,a,'j) is chosen in this form in order to represent (in the 
discretized version) the basic model of the one-component excitable medium, 
which is described by the Kolmogorov-Petrovsky-Piskunov type equation: 

dx{z,t) _ d^x{z,t) 

dt dz^ ^ ^ ' ■ 

This equation is frequently used in biophysical problems, in the combustion 
theory, in chemical kinetics (e.g. in the theory of the Belousov-Zhabotinsky 
reaction), solid-state physics, etc. (see, for example, [23,24,25,26]). One can 
see that the piecewise linear function G{x, a, 7) has a constant slope (1 — 2a). 
Depending on a the map x 1— > G{x,a,'~f) may exhibit qualitatively different 
behaviour. For < a < 1 (see Fig. la) it gives rise to a regular dynamics with 
one or two stable point, depending on the value of the parameter 7 (points A 
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and B in Fig. la), which are attractive for almost all points of the phase space. 
For a> 1 the slope of G{x,a,'~f) is larger then one (see Fig. lb). That means, 
that for the finite motion of phase points the dynamics would be chaotic. 

The uniform model (1) based on the coupled maps, which are described by 
functions (2), is known (see [22]). Following [22] consider this system as an 
iV-dimensional piecewise-linear map / : — > R^. Because G{x,a,'j) is 
the piecewise-linear function with the constant derivative, the differential D/ 
of the map / is a matrix with constant coefficients. This easily allows us to 
find the Lyapunov exponents Aj of such map. If ps are the eigenvalues of D/, 
then Xs = In \ ps\. If among all the values ps there exists an eigenvalue located 
outside the unit circle on the complex plane, then any trajectory of the map 
/ is unstable. Otherwise the model has the regular dynamics. Thus, to know 
qualitatively the behavior of the system it is sufficient to find the location of 
the eigenvalues ps on the complex plane. 

To find Ps we should analyse the characteristic polynomial for the differential 
of the A'^-dimensional map, which may be expressed via the determinants of 
three-diagonal matrices of different dimension. For these determinants recur- 
rent relations may be obtained. Taking the recurrent relations as difference 
equations with the given initial conditions we can find the solutions which 
are proportional to the Chebyshev polynomials with a linear function of the 
eigenvalues of D ^ [22] . Using then the properties of the Chebyshev polynomi- 
als it can be shown that the characteristic polynomial may be reduced to the 
quadratic trinomial of a certain function pg. Calculation of Ps and analysis of 
their location on the complex plane shows that the parametric space of the 
system (1), (2) has two regions with qualitatively different behavior (Fig.2). 

(1) The region Di, with the values of the parameters a and £, which cor- 
respond to the regular dynamics of the model. This implies that the 
modules of all the eigenvalues of the differential of the map / are less 
than unity. 

(2) The region D2 where at least one of the eigenvalues Ps satisfies the con- 
dition Ip^I > 1. Generally, the system may have an infinite dynamics. 
However, if the dynamics is finite the system (1), (2) behaves chaotically 
[9,10,11]. 

The system (1), (2) is an uniform since all its elements are identical (one has 
the same map x ^ G{x,a,'-f) for all the elements). Futher we will assume 
that each element of the system has the same map Xi 1— > G{xi, ai, 7) but along 
the chain the values of the control parameter ai may differ ^ . We assume that 
only two types of elements exist in the system. The corresponding values we 
denote as a and /3. As for the case of uniform system we will use the periodic 

^ In general, we can consider the case of different 7. But we will see that the 
dynamics of both uniform and nonuniform chains does not depends on 7. 
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boundary conditions, Xn+N{k)=Xn{k). 



2.2 The model with spatially periodic nonhomegeneity 



The periodic spatial nonhomegeneity imphes in our case that the chain has ele- 
ments with alternating values of the parameters in the map G, e.g. a, /9, a, /3, . . . . 
This particular nonuniform model allows a complete analytical analysis. 

Consider the system of coupled elements: 

Xn{k + 1) = G(x„(A;),Q!,7) + e{xnJ^i{k) - 2xn{k) + a;„_i(A;)), n is odd, 

< 

Xn{k + 1) = G{xn{k),l3,-f) + e(xn+i{k) - 2xn{k) + Xn-i{k)), n is even. 



where Xn+N{k) = Xn{k), N is even (to satisfy the periodic boundary condi- 
tion), and the function G is given by (2). To analyse the dynamics of this 
system we calculate the Lyapunov exponents using the approach described 
above. 

The differential of the A^-dimensional map corresponding to the system (3) 
reads 



N 



V 



2e-2a e 
s 1- 2s -2(3 
e 

e 





e 

1 - 2£ - 2a 







e 





l-2e-2(3j 



To find the eigenvalues p^, s = 1, 2, . . . , A'" of the matrix Q^, we calculate the 
determinant 



det{QN - pIn) 



2zie e ... £ 
£ 2z2e s . . . 
e 2zie e ... 

... £ 
e 0... 2z2e 
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with zi and Z2 defined by the relations l — 2e~2a — p — 2zie, 1 — 2s — 2/3 — p 
2z2S. Making expansion of the determinant we obtain: 



det{QN - rIn) = 2zieBN-i - e Bn-2 - 2s - s Bn-2, 



where 



Bm — 



2zie e 
e 2z2e e 
e 2zie 



. . . 2z2e 









Bat — 



2z2e e 
e 2zie e 
e 2z2e 



. . . 2zie 









It is not hard to see that Bm{z\, Z2) may be obtained from Biq{z\, Z2) exchang- 
ing the arguments zi ^ Z2. This suggest the way of finding Bn{zi,Z2) using 
the recurrent relation. Expanding Sjv (along the first line) we obtain: 



Bn = 2zisBn_i — s^Bn_2, N is even; 



= 2z2eBN-i — £^-Bjv-2, N is odd. 



This allows us to express determinants of the odd order via determinants of 
the even one, 



B. 



2N+1 — 



B2N+2 + £^B2N 

2z^ 



(4) 



and obtain then the recurrent relation for the even determinants 

B2N = S^(AziZ2 - 2)B2N-2 - £^B2N-4 ■ (5) 

One can treat Eq.(5) as a difference equation with the initial conditions 

Bo = 1, ^2 = £'(4ziZ2 - 1) ■ (6) 



Solving then (5), (6) together with (4), we arrive at: 

Bn{zi, Z2) — e^UM{z), for even N; 



Bn{zu Z2) = e^J^UNiz), for odd N, 
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where Un{z) is the second order Chebyshev polynomial. In turn, for Bf^{zi, Z2) 
we get: 

' Bn{zi, Z2) — £^U]s[{z), for even N ; 



Bn{zi,Z2) =e''J-UN{z), for odd N. 

Now, using the obtained relations for B^-i, -BAr-2 and we finally arrive 

at: 

det(g^ - pIn) = 2e''(zUN-i{z) - Un-2{z) - l). (7) 



This allows us to find the eigenvalues of the differential of the map (3), (2). 
As follows from Eq.(7), 

zUn-i{z)-Un-2{z)-1 = Q, (8) 



or, Tn{z) — 1 = 0. Here wc used the relation between the second order Cheby- 
shev polynomials Un{z) and the first order Tn{z). Since 



2 

where z'^ = Z1Z2, then with t = z — \J z^ — \ we obtain from Eq.(8): 

1 /.A, 1 



or (t^ - 1)^ = 0. Therefore, t, = e'^^, s = 1,2,...,N, and Zs = cos{2tts/N), s 

1,2, . . . , N. Taking into account that 2zie ■ 2z2£ = 'iz^e^, we obtain a simple 
equation: 

(1 - 2£ - 2q; - p)(l - 2£ - 2/3 - p) = Ae^ cos^ (jj'^ ' 

Solving it with respect to p, we find the set of eigenvalues of the differential 
of the map (3), (2): 



pl^2 ^l-2e-a-P±J{a-pr + 4e^ cos^ (^s) , (9) 



where s = 1,2,..., N/2. The dynamics of the system (3), (2) would be com- 
pletely regular if all the values of Ps were within the unit circle on the complex 
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plane. Since all the pf 2 real, the condition of the regular dynamics is: 



'2e + a + P> ^J4e^ + {a - (3)^, 



< 




< 1. 



The first inequality is always realized when s,a,(3 > 0. Thus, second one 
defines the restriction for the region with regular dynamics. 

In its turn, in the region of the parametric space where the Lyapunov exponent 



As = In > the condition e + {a + P)/2 + ^Je'^ + {a - /3)^/4 > 1 is satisfied. 
Therefore, the boundary separating two regions with regular and chaotic dy- 
namics is a surface in the three-dimensional parametric space {a, (3, e), which 
is given by the relation: 



To illustrate the obtained results wc show in Fig. 3 cross-sections of this surface 
by the planes 6 = const, with S = f3 — a being the nonhomegeneity parameter. 
Denote Di the region with regular dynamics and D2 the region where As > 
0. These two regions are separated by the boundary, defined (implicitly) by 
Eq.(lO). In Fig. 3 we also show by dashed line the corresponding boundary 
between Di and D2 for the case of homogeneous {5 — 0) system of coupled 
maps (Fig. 2). 

As it follows from the Fig. 3, for 5 > 1 for all values of e, a the model demon- 
strates the chaotic dynamics. For 5 < 1 the region Di arises where the system 
has regular dynamics. This region where < 5 < 1 is, in fact, a subregion 
of the corresponding region of the regular dynamics of the homogeneous 
model. 

For —1<S<0 there exist such e and a that the homogeneous model exhibits 
the chaotic dynamics, while the inhomogeneous still has the regular one. Note 
that the restriction a > —5 corresponds to the condition j3 > 0. For d < —1 
the region Di vanishes and for all e, a the model (3), (2) demonstrates only 
the chaotic dynamics. As it also follows from Fig. 3, two qualitatively different 
regimes make take place: a belongs to the region of the regular behavior, /3 
to the region of chaos with A^ > 0, with the global dynamics being regular. 
Similarly one can observe just the opposite case of the chaotic dynamics of 
the periodically inhomogeneous model. Realisation of the particular regime 
depends on the value of the diffusion constant e. 





(10) 
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2.3 Annular model with a single defect 



Let us now turn to the case of spatially inhomogeneous model of diffusively 
coupled maps for the case of a single "defect" . In this case among N element 
of the system only one element has the parameter /3 (see (2)), while the other 
(A^ — 1) elements have parameter a. Assume for simplicity that this element 
is located at n = 1. Then we can write for this model: 

Xn{k + 1) = G(xn(k),a,-f) +s(xn-i(k) - 2xn{k) +Xn+i(k)y n = 2, 3, . 

xi{k + l) = G{xi{k), + e(xN{k) - 2xi{k) + X2{k)), 



where Xn+N = Xn{k),n — 1,2, . . . , N. Let the function G{x, a, 7) is still given 
by (2). Calculate now the differential of the corresponding map for the system 
(11), (2): 



Df=QN^ 



l-2e-2p 
e 1 




e 

2e-2a 





e 

2e- 



2a ... 



e 



1 - 2£ - 2a 



.(12) 



If we apply the same approach as previously, to find the Lyapunov exponents 
of the model (11), (2), we should have to get the roots of the polynomial of 
the 2A^ + 2-th order. This can be done only numerically. Instead we will use 
somewhat different approach, which significantly facilitates the analysis. 

Applying the estimate of the eigenvalues of the matrix (12) according to the 
Gershgorin's theorem [27], we may evaluate the location of the region of the 
regular dynamics. By this theorem, all the eigenvalues of the matrix A — 
{%}nxn which belongs to the union of the discs on the complex plane: 

n 

Ps e \j{z e C : \z-aij\ < i?-} , 

i=l 

where s = 1,2, ... ,N and R'^ is either a row almost-norm of the matrix A, 
R'i = J2%^ Wijli OT a column almost-norm of the matrix A, R'- = X]'?=i Wijl- 

Since the differential of the map (12) is a real symmetric matrix then, (i) the 
row and column almost-norms are equal, and (ii) the eigenvalues of the matrix 
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D / are real. Therefore, the Gershgorin's discs transform on the real axis as: 

n 

Pse \j{zeR: \z-aij\ <R'i}, s^l,2,...,N . 



=1 

Because £ > 0, from (12) we get: 



n 



■ = ^ \aij\ = \e\ + \e\ = 2e . 



j=i 



Moreover, since docs not depend on i and in the differential of map Dj (12) 
there exist only two different diagonal elements, the eigenvalues of D j belong 
to the union of the two intervals: 

e R : |1 - 2£ - 2a - ^1 < 2£} \^\z e R : |1 - 2£ - 2/3 - ^| < 2£}. 

Then, taking into account that a, (5, e > 0, we obtain the upper hovLa<\ for the 
modules of the eigenvalues of Dj: 



< 



/3 + 2s, a< 13] 
a-\-2e, a > p. 



Thus, the region Di in the parametric space (a,/3, £), which satisfies the con- 
ditions 

r/3 + 2£<l, iia<(3] 
D, : (13) 
I q; + 2£ < 1, if q; > /3; 

corresponds to the regular dynamics of the inhomogeneous model (11), (2), 
while the region 

((3 + 2£>l, iia<^; 

D2 ■ (14) 
I q; + 2£ > 1, if a > 

corresponds to the positive A^. 

The regions Di and D2 are shown in Fig.4. The region Di is the lower bound 
of the region with the regular dynamics, because it is obtained via the upper 
bound of the eigenvalues of the map differential. This means that any point 
of Di would correspond to the regular dynamics, while any point of D2 would 
correspond either the regular dynamics, or (if the motion is finite) the chaotic 
one. 
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It is significant that the obtained estimate is apphcable not only for the model 
with a single defect, but for a wide class of annular chains of maps with 
diffusive coupling (1), which have only two types of elements (^with different 

parameters in the function ^(xj, «j, 7)^ . This follows from the fact that 
the particular number of elements of each type and their mutual location 
is not important. Indeed, almost-norms of the differentials of the maps of 
all such systems are the same and the centres of the two possible intervals, 
which contain according to the Gershgorin's theorem, the eigenvalues of Dj, 
do coincide. 



In the case when the two intervals, which we use to estimate ps, do not 
intersect, there is, according to the Gershgorin's theorem, "clusterization" 
of the eigenvalues of the differential of the map. Namely, q eigenvalues ps, 
s = 1,2, ... ,q, will belong to one interval, and {N — q) eigenvalues ps, s ~ 
q+l,q + 2, . . . , N will belong to the other interval. Here q is the number of el- 
ements with parameter a, and, respectively, {N — q) is the number of elements 
with parameter /?. Unfortunately, this does not help to improve the estimate 
of the region with the regular dynamics. 

Let us compare the estimate (13) with the exact result, that has been obtained 
above for the model with the periodic nonhomegeneity. Obviously, the latter 
model belongs to the class of models for which the this estimate holds true. 
To illustrate this we show the cross-sections 6 = const of the region Di, which 
corresponds to the regular dynamics of (3), (2), and of its estimate, the region 
Di (see Fig. 5). Prom Fig. 5 it follows that the estimate (13) of the region with 
the regular dynamics is rather adequate, however it does not show that for 
—1<S<0 the region of the regular dynamics of (3), (2) includes the region 
of the regular dynamics of the spatially homogeneous system. 

Turn now to the map (11), (2), which refer to the case of a single defect and 
find the eigenvalues of the differential of map (12) numerically. One can write 
for the characteristic polynomial: 



det{QN - pIn) 



2z£ + 2{a-(3) e ... e 
e 2ze e ... 

e 2ze ... 







2ze 



where 2ze — 1 — 2e — 2a — p. Expanding det{QN — pIn) in the first line we 
arrive at 



det{QN - pIn) = [2ze + 2{a - P)]Bn-i - 2e'Bj^_2 - 2(-l)^£^, (15) 
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with 



Bat — 



2ze e 
e 2ze e 
e 2ze 









2ze 



As previously, wc use the first fine factorization of -Bat to obtain the recurrent 
relation, Bj^ = 2zeB]\f_i — e'^Bjsi-2- Using the initial conditions, Bi = 2zs, 
B2 = £^(4z^ — 1), we obtain: 



Bn^e^Uniz), 



(16) 



where Un{z) is the second order Chebyshev polynomial. Taking into account 
Eq.(16) and representation of Un{z), 



Un{z) 



/^x iV+l 

+ v/i^^)^+i -{z- v/i^^)^+i _\-t) - ^""^^ 



1 

--t 
t 



where t — z — y/ z'^ — 1, we recast Eq.(15) into the following form: 



det{QN - pIn) = 2e 



N 



T-(-l) 



N 



2t^ \ 6 J {1- t^)t^ 

Thus, the characteristic equation for the differential of the map (12) reads 

' a — P 



,27V+2^0 (^Jl\ ^2iV+l_^2iV_2(_i)iV^iV+2^2(-l)^i^+i2-2 



t-1 - 0. 



Finding numerically the roots of this equation, wc obtain ps- Comparing ps 
with unity for different sets of a, /3, e, for the annular chain of diffusively 
coupled maps with a single defect we find the regions Di and D2 with qualita- 
tively different dynamics. In Fig. 6 we show cross-sections of these regions by 
the planes with the fixed values of 5 = const, where 6 = (3 — a. The estimate 
(13) of the region with the regular dynamics is marked by the dashed line. As 
in the case of the periodic nonhomegeneity, the actual region for positive 6 is 
larger than that predicted by the estimate. However for the negative 6 the es- 
timate is in a perfect agreement with the numerical results as well as with the 
results corresponding to the homogeneous chain. Thus, if one of the elements 
of the chain had the parameter smaller than that of the others, there would 
be no impact on the global dynamics of the system. On the other hand, if the 
parameter which refers to the defect had larger value than that of the other 
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elements, this would diminish the region of the regular dynamics as compared 
to the case of the uniform system. 



3 Conclusions 



We investigated the dynamics of diffusively coupled one-dimensional maps, 
which model one-dimcnsional excitable medium with space nonhomegeneity. 
We considered different types of the nonhomogcncitics. For all systems, with 
the spatial nonhomegeneity, which may be described in terms of two types 
of elements with different values of the control parameter, we found analyti- 
cally a lower bound of the region in the parametric space where the ensemble 
exhibits the regular dynamics. This estimate has been obtained by localizing 
the eigenvalues of the differential of the map. For two different model of this 
class we also obtained exact analytical and numerical results which completely 
characterize the dynamics of the system. 

We observed that all three approaches (the approximate analytical analysis of 
the chain with two types of elements (13), the exact solution for the spatially 
periodic defects (§2.2) and numerical study (§2.3) of the model with a single 
defect) show the possibility of two qualitatively different types of the ensemble 
behavior. In the first case the system demonstrates a regular global behavior 
when one value of its parameter belongs to the region with the regular dy- 
namics and the other one belongs to the region with the unstable dynamics. 
Vice versa, in the second case the global dynamics occurs to be unstable, and 
if the motion is finite — chaotic. The particular type of the chain behavior de- 
pends on the diffusive parameter e. Obviously, if £ is large enough, the global 
dynamics would not be regular. 

We wish to stress that according to our knowledge, the present study is the 
first attempt of the analytical description of the dynamics of inhomogeneous 
systems of coupled maps. These systems, being of significant importance from 
the point of view of applications, may exhibit a rich physical behavior, giving 
rise to a lot of new and interesting problems. One can mention such prob- 
lems as properties of the wave travelling along the chain and influence of 
the defects on these motions, dynamics of the ensemble with time-dependent 
parameters, controllability conditions, synchronization and many others (see, 
e.g., [29,30,31,18,19,20]). These problems, not discussed here will be addressed 
elsewhere [28,32]. 
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a) b) 

Fig. 3. Cross-sections by the planes 5 = const of the parametric space of the model 
(2), (3) for -1< (5 < (a) and < 5 < 1 (b). 
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a) b) 

Fig. 5. The estimate and the exact result for the region with the regular dynamics 
for the model with periodic nonhomegeneity for — 1 < 5 < (a) and < 6 < 1 (b). 
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Fig. 6. Cross-sections by the planes 6 = const of the parametric space for the chain 
of = 20 elements with a single defect for 5 = —0.3 (a) and 5 = 0.3 (b). 
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